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Abstract. In the inflationary paradigm it is believed that quantum fluctuations during inflation 
seeded the temperature anisotropies seen in the CMB. We present an extension of the semiclas- 
sical Einstein equations which couples n-point correlation functions of a stochastic Einstein 
tensor to the n-point functions of the quantum stress-energy tensor. To apply this extension to 
calculate the quantum fluctuations during an inflationary period, we take as a model a massive 
conformally coupled scalar field on a perturbed de Sitter space and describe how a renormaliza- 
tion independent, almost-scale-invariant power spectrum is produced. Furthermore, we discuss 
how this model yields a natural basis for the calculation of non-Gaussianities of the curvature 
fluctuations. 

1. Introduction 

Shortly after the discovery of the cosmic microwave background (CMB) by Penzias and Wilson 
[32], Sachs and Wolfe predicted anisotropies in the angular temperature distribution [37]. In 
their famous paper they discuss what was later coined the Sachs-Wolfe effect: The redshift in 
the microwave radiation caused by fluctuations in the gravitational field and the corresponding 
matter density fluctuations. As a next step it is prudent to ask how the fluctuations imprinted 
upon the CMB were created, i.e., to explain the origin and predict the form of the initial curvature 
fluctuations. In the inflationary paradigm it is claimed that the primordial fluctuations are seeded 
by quantum fluctuations during inflation [14, 16, 28, 29]. 

The usual computation of the power spectrum of the initial fluctuations produced by single- 
field inflation can be sketched as follows [6, 14, 16]: First, one introduces a (perturbed) classical 
scalar field tp + dtp, the inflaton field, which is coupled to a (perturbed) expanding spacetime 
g + Sg. Then, taking the Einstein equation and the Klein-Gordon equation at first order in the 
perturbation variables, one constructs an equation of motion for the Mukhanov-Sasaki variable 
Q — dtp + tpH l <t>, where <E> is the Bardeen potential [5] and H the Hubble constant. Q is 
then quantized 1 (in the slow-roll approximation) and one chooses as the state of the associated 
quantum field a Bunch-Davies-like state. Last, one evaluates the power spectrum Pg(k) of Q, 
i.e., the Fourier-transformed two-point distribution of the quantum state, in the super-Hubble 
regime k <K aH and obtains an expression of the form 2 



where Aq is the amplitude of the fluctuations, ko a pivot scale and n s the spectral index. Notice 
the factor of AT 3 in (1.1) which gives the spectrum the 'scale-invariant' Harrison-ZeV dovich 
form if n s — 1 . Depending on the details of model, n s < 1 and there is also a possibility for a 
scale dependence of n s - the 'running' of the spectral index n s = n s {k). 

This result can then be related to the power spectrum of the comoving curvature perturbation 
% which is proportional to Q and compared with observational data. Assuming adiabatic 
and Gaussian initial perturbations, the WMAP collaboration finds n s = 0.9608 + 0.0080 (at 
ko = 0.002 Mpc~') in a model without running spectral index and gravitational waves, excluding 

1 A recent discussion about the quantization of a such system can be found in [17]. 
2 An alternative definition of the power spectrum is PqQc) = (2n)~ 2 k 3 Pq(U). 
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a scale-invariant spectrum at 5<x [20]. Furthermore, the data of WMAP and other experiments 
can be used to constrain the deviations from a pure Gaussian spectrum, the so called non- 
Gaussianities, that arise in some inflationary models [6, 7, 25]. 

In [1, 2, 3, 31] concerns have been raised whether the calculation leading to (1.1) and similar 
calculations are correct: The authors argue that the two-point distribution of the curvature 
fluctuations has to be regularized and renormalized similarly to what is done in semiclassical 
gravity. As a result the power spectrum is changed sufficiently that previously observationally 
exluded inflation models become realistic again. On the contrary the authors of [15, 26] argue 
that the adiabatic regularization employed in [1, 2, 3, 31] is not appropriate for low momentum 
modes if evaluated at the Horizon crossing and irrelevant for these modes if evaluated at the end 
of inflation. We will come back to the issue of renormalization later in this paper and find that in 
our model the curvature fluctuations must not be renormalized. 

A slightly different approach to the calculation of the power spectrum based on stochastic 
gravity can be found in [23, 35, 36]. In spirit similar to the approach presented in this paper, 
the authors equate fluctuations of the stress-energy tensor with the correlation function of the 
Bardeen potential. In the super-Hubble regime they obtain an almost scale-invariant power 
spectrum. Moreover, they discuss the equivalence of their stochastic gravity approach with the 
usual approach of quantizing metric perturbations. 

Inspired by these works we follow an approach strictly different from the standard one. 
Instead of quantizing a coupled system of linear inflaton and gravitational perturbations, we 
aim at extending the semiclassical Einstein equation to describe curvature fluctuations via the 
fluctuations in the stress-energy tensor of a quantum field. More precisely, we choose as quantum 
matter a scalar quantum field and analyze carefully the expectation values of the corresponding 
stress-energy tensor and of their products over a Gaussian Hadamard state. We then interpret the 
Einstein tensor as a stochastic field and equate its n-point function with the symmetrized «-point 
function of the quantum stress-energy tensor. Subsequently we use the obtained system to discuss 
the influence of quantum fluctuations on the Bardeen potential over the flat patch of a de Sitter 
background. We prove that in the limit of early time the power spectrum of these perturbations 
tends to the Harrison-Zel'dovich spectrum. The results show that no special features of the 
metric nor any renormalization ambiguity enters the computation; only the Hadamard form of 
the underlying state matters. Finally we demonstrate that, despite the linearity of the quantum 
field, the three-point function of the Bardeen potential does not vanish. 

The paper is organized as follows. In the next section we shall briefly discuss the form 
of the matter quantum field we shall use and its backreaction on the background by means of 
semiclassical Einstein equation. The basic idea about the influence of quantum matter fluctuations 
on curvature perturbations is outlined in the third section. The fourth section contains the 
construction of a simple model on de Sitter spacetime, the analysis of the power spectrum of 
certain fluctuations and a brief discussion on the naturally obtained non-Gaussianities. Finally, 
the last section, contains a discussion and a summary of the obtained results. 

2. Review: Semiclassical Einstein equation and scalar fields 

The semiclassical Einstein equation 3 

G ab = u(:T ab :) (2.1) 

for a quantum state u> and a (conserved) quantum stress-energy tensor T a i, on a globally hyperbolic 
spacetime (M, g) allows us to analyze the backreaction of quantum fields on curvature, see e.g. 
[11,40,41]. 

In this paper, for simplicity, we will take as the quantum field a free conformally coupled 
scalar field ip with mass to, i.e., the classical dynamics are governed by the Klein-Gordon 

3 We use units where c = ft = 8jtG=l. 
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equation Pip = —Dip + lR<p + m 2 ip - 0. The stress-energy tensor associated to this scalar field 
can be written as 

2 1 1 / 2 R \ 

Tab = ^(d a (p)(db<p) - -gab(d c (f)(d c tp) - -<p \m g ab + -g a b - R a b + 2V a d b j <p, 

which differs from the standard definition by the term \g a b f> P<P- This term was introduced in 
[27] to account for the trace anomaly in the quantum theory. A quick calculation gives us the 
trace of T a b as 

T = tf b T ab = -n?<? + -tpVq>. (2.2) 

The quantization of scalar fields on globally hyperbolic manifolds in the algebraic approach 
has been discussed thorougly in the existing literature, see e.g. [9], and will not be repeated here. 
To accomplish the quantization we shall assume that the quantum state to is a quasi-free state, 4 
hence all its zi -point distributions oj„ descends from its two-point distribution 0)2. Furthermore, 
to make (2.1) well-defined, we will require that a>2 satisfies the microfocal spectrum condition 
which is equivalent to it being of Hadamard form [8, 34, 38], i.e., if x, y e M are in a geodesically 
convex neighbourhood, a>2(x, y) is of the form 

to 2 = lim I— +Vln — | + W = ( H + W, (2.3) 

s^\cr s A) 

where cr s (x, y) = a(x, y) - is (t(x) - t(y)) + e 2 with <x being the squared signed geodesic distance 
and t a time function, U(x, y) and V(x, y) = V„(x, y) cr(x, yf are smooth functions depending 
only on the local geometry and W(x, y) is a smooth, symmetric function which characterizes the 
state. 

The distribution H in (2.3) is called the Hadamard singularity and leads to the singular UV 
behaviour in the quantum theory. When we regularize a>, by point-splitting regularization, to yield 
a finite Einstein tensor on the left-hand side of (2.1), 'H is subtracted from cl>2- This regularization 
is indicated by the normal ordering : ■ : in (2.1), i.e., in the semiclassical Einstein equation we are 
evaluating Wick polynomials. The ambiguity in the normal ordering prescription leading to the 
renormalization freedom has been rigorously analyzed in [21, 22], where it was shown that in a 
locally covariant quantum field theory the ambiguity amounts to adding polynomials in local 
curvature and the mass parameter. 

We thus note the following regarding the trace of the semiclassical Einstein equation (2.1): 
Selecting a Hadamard state to for the scalar field <p on the globally hyperbolic spacetime (M, g), 
we obtain (see also [11, 21]) 

-R = to(:T:) = -m 2 [W] + 2[Vi] + am 4 +f5m 2 R + y OR, 

where [•] denotes the Synge bracket [f](x) = fix, x), [W] is a state dependent contribution, [Vi] 
is due to the trace anomaly [42], and a,f3,y are renormalization constants. Since we want to 
have a well-posed initial value problem, according to [41], we shall choose y — 0. Furthermore, 
different choices of /3 in the equation —R - oj(:T:), can be reabsorbed in the redefinition of the 
Newton constant. Because we do not want to change it, we shall fix [5 — 0. With these choices 
we notice that the single remaining freedom is in the constant a whose value also depends on the 
choice of the scale A in the Hadamard singularity (2.3). 

3. Higher moments of the semiclassical Einstein equation 

As noted in the introduction above, in the semiclassical approach we are equating a classical 
quantity, the Einstein tensor, with the expectation value of a quantum observable, the quantum 
stress-energy tensor, i.e., a probabilistic quantity. Such a system could make sense only when 

4 A quasi-free state is also called Gaussian state. 
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the fluctuations of the quantum stress-energy tensor can be neglected. Unfortunately, as also 
noticed in [33], the variance of quantum unsmeared stress-energy tensor is always divergent 
even when proper regularization methods are considered. The situation is slightly better when a 
smeared stress-energy tensor is analyzed. In that way, however, the covariance of (2.1) gets lost. 
A possible way out is to allow for random perturbations also on the left hand side of (2.1). This 
is the point of view we shall assume in this paper. Notice that this is similar in spirit with what 
happens in the analysis of the random forces leading to Brownian motions by means of Langevin 
equations. 

Consider now the Einstein tensor as a random field. Then we could imagine to equate the 
probability distribution of the Einstein tensor with the probability distribution of the stress-energy 
tensor. This suggestion, however, seems largely void without a possibility of actually computing 
the probability distributions of the stress-energy tensor because, as discussed above, its moments 
of order larger than one are divergent. 

Instead we may approach this idea by equating the hierarchy of n-point functions of the 
Einstein tensor with that of the stress-energy tensor: 

{G ab {xi)) = to(;.T ab {x\y), (3.1a) 
(6G M 6G C , AX 2)) = \ <*:6TM): :c^(* 2 ): + UST„M): :6T ab ( Xl ):), (3.1b) 

((6Gf n (x 1 ,...,x„)) = &[a>(:6T:® n (xi,...,x n ))], n>\, (3.1c) 

where a> is a Hadamard state and we defined 6G a b = G a b - (Gab) an d '-ST a b'- - '-Tab'- - w(:T a i,:). 5 
The symmetrization on the right-hand side, which we denote by (5, is necessary because the 
classical quantity on left-hand side is invariant under permutation of the points x, . 

We emphasize that we are equating singular functions in (3.1c). Having all the n-point 
functions of the Einstein stochastic tensor, we can easily construct an equation for the moments 
of the smeared Einstein tensor which equals the moments of a smeared stress-energy tensor by 
smearing both sides of (3.1) with tensor products of a smooth compactly supported function. 
This smearing also automatically accounts for the symmetrization in (3.1). 

Furthermore we stress that equating moments, obtained smearing both side of (3.1), is not 
equivalent to equating probability distributions. Although it is also possible to arrive at a de- 
scription in terms of moments when coming from a probability distribution, the inverse mapping 
is not necessarily well-defined. Successful attempts to construct a probability distribution for 
smeared stress-energy tensors can be found in [18, 19]. 

Consider now a Gaussian Hadamard state to of a conformally coupled scalar field (sona 
spacetime (M,g), the background spacetime. Our aim is to calculate the perturbation of the 
background spacetime as specified by the correlation functions on the left-hand side of (3.1) due 
to the fluctuations of the stress-energy in the quantum state to as specified on the right-hand side 
of (3.1). In particular we will require that u satisfies (3.1a) when we identify the Einstein tensor 
of the background spacetime G„b with (G a b). Note that by choosing this ansatz we are completely 
ignoring any backreaction effects to the background metric and evaluate the stress-energy tensor 
on a state specified on the background spacetime. 

Later on we consider curvature perturbations induced by a simple metric perturbation. For 
this reason it will be sufficient to work with the trace of (3.1) (using the background metric) 
instead of the full equations. With the definition 

ci —ah,-* 

o = -g U a b 



5 The symbol Kl stands for the tensor product at different points, the so called exterior tensor product. 
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Figure 1. A few graphs illustrating (3.2c) for n — 2, n — 3 and n — 5. 



the equations (3.1) simplify to 



(S(xO) = m 2 [W] - 2[Vi] - am 4 , 
(S (jfi) 5 (x 2 )> - <S (xi))<5 (jc 2 )> = m 4 (w2(JCi» x 2 ) + io\(x 2 , x x )), 



(3.2a) 
(3.2b) 



<(S -{S)T n {xi x n )) = 2 n m M e 



Zi— r U) 2 (Xi, Xj) « 
1 1 jC~\ 



n > 1, 



(3.2c) 



where the sum is over all directed graphs G with « vertices 1 , . . . , n with two arrows at every 
vertex directed to a vertex with a larger label. Af. e {0, 1,2} is the number of arrows from i to 
j. If we perform the symmetrization in (3.2c), we see that the sum is over all acyclical directed 
graphs with two arrows at every vertex. For illustration some graphs are shown in Fig. 1. 

To obtain (3.2b) and (3.2c), note that :ip 2 : - u>(:ip 2 :) doesn't depend on the choice of normal 
ordering 6 and thus only (3.2a) needs to be renormalized. Therefore we may choose normal 
ordering with respect to <jj 2 t° see that the combinatorics are equivalent to those in Minkowski 
space. Moreover, as u 2 is a bisolution of the Klein-Gordon equation, the term \ fP(fi which 
causes the conformal anomaly in (3.2a) does not contribute to the higher moments. 



4. Fluctuations around a de Sitter spacetime 

We shall now specialize the general discussion presented above to a special class of fluctuations 
around an exponentially expanding, flat FLRW universe. That is, the background spacetime 
(M, g) is given in conformal time t < by the metric tensor 

g = (flVT 2 (-dr ® dr + <5y dx l ® dx J ) 

and we consider fluctuations derived from metric perturbations of the form 

g = (Ht)- 2 ( - (1 + 2Y) dr ® dr + (1 - 2»F) 8 tj Ax 1 ® Ax'). (4. 1) 

The kind of fluctuations that we consider by choosing (4. 1) resemble those that are present in 
single-scalar field inflation in the longitudinal gauge, where there are only 'scalar fluctuations' 

6 Indeed this holds true if we replace if 1 with Up 1 , for any linear operator L. 
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without anisotropic stress (so that the two Bardeen potentials coincide = [16, 29]. Notice 
that, for classical metric perturbation, these constraints descend from the linearized Einstein 
equation, however, a priori there is no similar constraint in (3.1b). Despite these facts, we 
proceed analyzing the influence of quantum matter on this special kind of metric perturbations 
and we also refrain from discussing the gauge problem associated to choosing a perturbed 
spactime. 

We can now calculate the various perturbed curvature tensors and obtain in particular 

S = 12H 2 (1 - 3W) + 24ff 2 r v P' - 6H 2 t 2x V" + 2H 2 t 2 A¥ + OCP 2 ) 

for the trace of the perturbed Einstein tensor, where A is the standard Laplace operator. Dropping 
terms of higher than linear order, this can also be written as 

S - (S) = -6H 2 t 4 ^ - ^aJt- 2 ^, (4.2) 

where (S) = VLB 2 is nothing but the scalar curvature of the background spacetime. Notice that, 
up to a rescaling, the operator on the right hand side of (4.2) looks like a wave operator with the 
characteristic velocity equal to 1 / V3 of the velocity of light. 

We can now evaluate the influence of quantum matter fluctuations on the metric fluctuations 
by inverting the previous hyperbolic operator by means of its retarded fundamental solutions A# 
and applying it on both sides of (3.2b) and (3.2c). From (3.2b) we can then obtain the two-point 
correlation functions of ¥ (per definition (W) = 0): 

(SP(jfi)¥(jC2)> = m 4 ff A R (xuyi)A R (x2,y2)(col(yi,y 2 ) + o)l(y2,yi))d 4 y 1 d' l y2. (4.3) 
JJr 8 

Employing the retarded fundamental solutions in the inversion without adding any solution of 
(4.2), we are implicitly assuming that all the «-point functions of the Bardeen potential *F are 
sourced by quantum fluctuations. Here we are only interested in evaluating their effect. 



4.1. Bunch-Davis state and the squared two-point distribution 

In order to proceed with our analysis, we shall specify the quantum state a> for the matter theory. 
Following the ansatz discussed in the preceding section, we choose a quasi-free Hadamard state 
which satisfies the semiclassical Einstein equation on the background. In particular, we require 
that a> solves (3.2a), namely 

12H 2 = m 2 [W] - 2[Vi] - am 4 . 

The right hand side of the previous equation is characterized by three contributions: The state 
dependent part m 2 [VK], the anomaly part 2[Vi], which takes a very simple form and is proportional 
to H 4 , and the renormalization freedom am 4 (remember that we fixed the other renormalization 
constants to zero). For the semiclassical Einstein equation to hold, we therefore have to require 
that [W] is a constant. Then, having fixed H and m (no matter their absolute value), there is 
always a choice of a for which the chosen metric g and a> satisfy the semiclassical Einstein 
equation. On a de Sitter spacetime these criteria are satisfied by the well known Bunch-Davis 
state cl>bd [10] - the only Hadamard state which is invariant under the symmetry group of de 
Sitter space. 

In order to evaluate the influence of the quantum matter fluctuations on the Bardeen potential 
via equation (4.3), we have to discuss the form of the two-point distribution of the chosen state 
and its square. The two-point distribution of the Bunch-Davis vacuum for a massive, conformally 
coupled field on de Sitter spacetime, takes the well known expression [4, 10, 24, 39] 



/3 3 „ l+Z\ /l 

ubd = ~z : r iF i - + v, - - v; 2; 



16jtcos(jtv) \2 '2 ' ' 2 ' V4 H 2 ' 
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where 2 F\ is the ordinary hypergeometric function and the necessary e-prescription has been 
omitted. Moreover, the auxiliary function Z{x\,Xz) appearing in the preceding formula is nothing 
but the geodesic distance in the five dimensional Minkwoski space into which de Sitter space can 
be embedded as an hyperboloid. Omitting again the necessary e-prescription, cobd can also be 
recast into 

H 2 , m 2 (3 3 1-Z\ 1-Z m 2 - 

— 8^ (1 - Z) + T6^'(2 +V '2- y;2; — + J^ W 

with a certain smooth function W(Z{xi , x 2 )). 

In the conformally flat patch of the de Sitter spacetime Z takes a simple form and can be 
expressed in terms of the squared, signed geodesic distance of the conformally related Minkowski 
space as 

Z(x u x 2 ) = 1 - °~ M ( Xl,Xz \ cr M Oi,x 2 ) = -(Ti -t 2 ) 2 + (x l - x 2 ) 2 . 
2tit 2 

On (M, g) we therefore obtain an expression for lobd which resembles the Hadmard form on the 
conformally related Minkowski space: 

<jj B d(x\,x 2 ) = lim — - — +less singular terms. (4.4) 

e->o+ 4ji z o~m(xi i x 2 ) + 2is (ri - t 2 ) + s L 

Notice that the less singular contributions vanish in the limit of zero mass. 

As can be seen in (3.2b), we need to compute the square of two-point distribution of the 
state in question, i.e., in this case the square of u>bd- For our purposes it will be sufficient to 
compute the square of the leading singularity in the Hadamard state. To this end we notice that, 
up to a trivial rescaling by H 2 t\t 2 , the leading singularity coincides with the vacuum state com 
for a massless field theory on the conformally related Minkowski spacetime. The square of the 
massless two-point distribution on Minkowski space is 

\2 



Um(xi,x 2 ) = lim 



1 1 



e->6+ \4n 2 ct m (x 17 x 2 ) + 2\e{t\ - t 2 ) + e 2 , 

Writing wm in terms of its spatial Fourier transform, an expression for the spatial Fourier 
transform of the square of the massless Minkowski vacuum can be obtained as 

00 

<«&(jci,j&)= lim — ^- f e*ca-*> f e-^-^e^dpdf, (4.5) 

K- 1 k 

where xt = (t,, x/) and k = \k\. Later on we will use this expression in order to obtain the power 
spectrum of the fluctuations of 4 P. 

4.2. Power spectrum of the metric perturbations 

We want to compute the power spectrum P(t, k) of the two-point correlation of T* at the time r. 
Since both the spacetime and the chosen state are invariant under spatial translation, it can be 
defined as 

<'P(T,f 1 )'P(T,f 2 )) = -^-3 f P(r,£)e*«-^d¥. 

R 3 

To obtain P, we first need an expression for the retarded operator Ar corresponding to (4.2): 

T 

(A* f)(j, x) = -^-3 j J A R (r, Tl J) /(n , k) e i£f dr' d% with 

R3 -00 

A R (T,T U k) = -—p;^— sin 



m 2 r\ k \ V3 
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where / is a compactly supported, smooth function, ■ is the spatial Fourier transform. 7 We can 
then rewrite the spatial convolutions present in (4.3) as products in Fourier space to obtain 



T T 

P(t, k) = m 4 J J A r (t, n , k) A r (t, t 2 , k) 0) 2 bd {t\ , t 2 , k) At x dr 2 . 



— oo — oo 



Note that the symmetrization of the state is taken care of indirectly by the equal limits of the two 
integrations. 

As discussed above (see (4.4) and the following paragraph), we will compute the contribution 
due to the leading singularity of the Hadamard state: 



Po(x,k) 



T T 

H 4 m 4 J j A R 



(r, n , k) A r (t, r 2 ,k) t\ t\ «L(ti , r 2 , k) dr x dr 2 . 



We emphasize at this point that, because of the form of (4.5) and of A R , no ^-infrared singularity 
appears in Po(t, k) at finite r. Recall also that the error we are committing, using Po(t, k) at 
the place of P(t, k), tends to vanish in the limit of small masses. Inserting the spectrum of to 2 ^ 
obtained in (4.5) and switching the order in which the integrals are taken (for s > 0), we can 
write 

2 



-> m 4 r 1 
P (r,k) = lim ^ - 



where we have introduced the auxiliary function 



-Mr,-,. 



<r sp dp, 



(4.6) 



C KT 2 

A(T,K,p)= -j- 



sm(K(T-Ti))e- ipTl dn, 



(4.7) 



which can also be written in closed form in terms of the generalized exponential integral 8 £ 2 as 
A(t, k,p)= 1 -kt (E 2 (i (p + k)t) e 1KT - E 2 (i (p - k)t) & -' iKT ) 

for p > k > and by the complex conjugate of this expression if k > p,K > 0. In the following 
study of the form of the power spectrum Pq the auxiliary function A will be instrumental. 

Lemma 4.1. For \p\ + k > 0, Air, k, p) has the T-uniform bound 



4 k 2 
\k 1 - p l \ 



For large negative times it satisfies the limit 



(4.8) 



lim \A\ = -— 



Proof. Using the fact that 



-l/JTi 



+ K~ 



-lpTl 



K 2 — p 2 ' 



7 The conventions for the Fourier transform observed here are: 



/(*)= I /We- M d 3 f, /W = (2; 



f(k) e ' k -'d 3 k. 



(4.9) 



For a definition and various properties of these special functions see e.g. [30, Chap. 8]. 
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we can perform two integrations by parts to obtain 



A(t, k, p) = =■ (e^ T + R(t, k, p)), with 

a: z — 

T 

/' 4 6 
-7 cos (k(t - n)) H 7 sin(/<-(T - n)) 
Ui J 



e-" ,Tl dr, 



It is now easy to obtain an upper bound for R which is uniform in conformal time, namely \R\ < 3, 
which then yields the bound (4.8). 

For the second part of the proposition we perform a change of the integration variable to 
x = n/r: 

OO 

R(t,k, p ) = - rf-jcos(*T(l-Jc))+ — —7 sin (* t (1 - *)) | e^™' dx. 

J \X J K T X J 

1 

The contribution proportional to 1/r in R is bounded by C(ac)/|t] and thus vanishes in the limit 
r — > -oo. Moreover, since \p\ + k and 1/x 3 is L 1 on [1, oo), we can apply the Riemann-Lebesgue 
lemma and see that this contribution vanishes in the limit r — > -oo. The remaining part of \A\ is 
k 2 \k 2 - p 2 | _1 , which is independent of r, and thus the limit (4.9) holds true. □ 

Note that the bound for A obtained above is not optimal. Numerical integration indicates that 
|A| 2 is monotonically decreasing in r and thus bounded by the limit stated in (4.9) (see also Fig. 
2). Nevertheless, we can use this lemma to derive the following bounds and limits for Pq. 

Proposition 4.2. The leading contribution Pq to the power spectrum of the Bardeen potential *F 
induced by a conformally coupled massive scalar field in the Bunch-Davis state is bounded by 
the Harrison-Zel'dovich spectrum uniformly in time, namely 

i , Ai 16 C „ 3 - 2 V3 arccoth V3 4 

\Po(t, k) \< C = — m\ 

1 1 |^|3 384jt 2 

and it tends to the Harrison-Zel'dovich spectrum for t — * — oo, i.e., 

C 

lim P (T,k) = —. 

T^-OO | fe | 3 

Proof. The proof can be easily obtained using the r-uniform estimate (4.8) obtained in lemma 
4. 1 and computing the integral 



oo 2 



3 - 2 V3 arccoth V3 m 4 
24jt2 W 



k 

Having shown the first part of the proposition, let us now analyze the limit 

2 



lim P (tJ) = -^-r [ lim a(t,^-, P 

k 



dp, 



where we have taken the r-limit before the integral and already evaluated the e-limit because \A\ 2 
is bounded by an integrable function uniformly in time. Inserting the limit (4.9) from lemma 4.1, 
we can compute the /^-integral 



oo 



thus concluding the proof. 



: 3 - 2 V3 arccoth V3 m 4 
32n 2 J \3p 2 -k 2 j ~ r ~ 384ji 2 ~k? 

k 
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We can complement the results of proposition 4.2 with the following observation: 



Proposition 4.3. The power spectrum Po has the form 

Po(T,k) = — 5 — , 

l*p 

where Pq is a function of\k\r only. 

Proof. Noting that A(r, k, p) is a function of r k and t p only and performing the e-limit in (4.6) 
inside the integral, this can be seen by the substitution x — p r in (4.6). □ 

We would like to improve the estimate of Po(j, k) for t close to zero. Adhering to our 
previous strategy, we shall first give a new estimate for A(t, k, p): 

Lemma 4.4. The auxiliary function A(t, k, p) is bounded by 

tk 2 

\A(t,k,p)\<-2—, p + Q, t<0. 
\p\ 

Proof. Recalling the form of A given in (4.7) and integrating by parts, where we use that 

e-^' = ip- 1 d Tl Q- ipT \ we find 



A(t, k, p) 



\k 2 t 2 r 
p J 



-r COS(*(t- Tl)) + rSin^CT-T,)) 

— oo 



e 



We then take the absolute value and estimate the trigonometric functions, which gives us a bound 
on A, namely 



2 " 2 nT-Tl^ 



dn = -2— . 



Performing the integration in p analogously to the second part of proposition (4.2), the last 
lemma immediately leads to a corresponding bound for P : 

Proposition 4.5. The leading contribution Pq of the power spectrum of the fluctuations of the 
Bardeen potential *P satisfies the inequality 

\Po(r,k)\ 



72jt 2 |^| 

and therefore, in particular, Pq(0, k) — 0. 

The rescaled power spectrum !Po(|£|t) can be analyzed numerically and a plot is shown in 
Fig. 2. It clearly exhibits the asymptotic behaviour of Pq discussed in propositions 4.2 and 4.5. 

In this section we have used the leading singularity 9 of the two-point function of the Bunch- 
Davis state on a de Sitter universe to compute the influence of quantum matter on the power 
spectrum of the fluctuations of Bardeen potential. We have seen that this results in an almost 
scale-invariant power spectrum. We stress that such a singularity is not a special feature of the 
Bunch-Davis state but is common for every Hadamard state. Moreover, although our analysis 
has been done on a de Sitter universe, similar quantum states have been constructed on universes 
which are asymptotically de Sitter spaces in the past [12, 13]. All these states tend to the 
Bunch-Davis state for t — > -oo and are of Hadamard form. 



'Recall that considering only the leading singularity in the Bunch-Davis state also corresponds to the limit of 
vanishing mass. 
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Figure 2. Plot of the rescaled power spectrum !Po(|£|i"), where C is the same proportionality 
constant as in proposition 4.2. 



4.3. Non-Gaussianities of the metric perturbations 

It follows from (3.2c) that the n-point correlation for the Bardeen potential *P will, in general, 
not vanish. Also for odd n they will be different from zero and hence *F is not a Gaussian random 
field. As a first measure of the non-Gaussianity of one usually calculates its three-point 
correlation function or the corresponding bispectrum B: 



CP(T,fl)1'(T,f 2 ) < P(T,f3)) = 



1 



(2ji) 9 



6(k\ +Ic2 + h) B(t, k\ , ka, £3) 

x e**-*-*-*- 1 *-^ d 3 li d¥ 2 d%. 



Assuming non-zero k\, k\ and £3, we will derive the form of the bispectrum B considering 
(as above) only the contribution due to the leading singularity of the Bunch-Davis state, which 
we will denote by Bq. We will follow the same steps that lead us to the calculation of the power 
spectrum in the previous section. That is, we apply the retarded propagator A« of (4.2) as in (4.3) 
to the right-hand side of (3.2c) for n — 3 to obtain an equation for and insert for the two-point 
distribution the conformally rescaled two-point distribution of the massless Minkowski vacuum. 
The result can again be expressed in terms of the auxiliary function A defined in (4.7): 



Bq(t, hi,ki,ks) — lim 



0+ ia^l>k\k\k\ 



I 



' & -e (w, 7 (-£,)+w,,(£ 3 )+|p|) 

ojpi-k^cojjiki)^ 



vr 



r, —^,o)p(-ki) + \p\ 



x A 



t, — ,-o)p(k 3 ) - \p\ 
V3 



\k 2 \ 

t, —,a)p(h) - top(-ki) 
V3 



+ permutations 



d% (4.10) 



where u>p(k) — \k + p\ and the sum is over all permutations of 1, 2, 3. 

We can apply the same bound on A which has been used in the previous section to bound the 
power spectrum Pq to produce a bound on the integrand of Bq almost everywhere. 10 Nevertheless, 



10 We cannot bound the integrand of Bo everywhere using (4.8) because \k2\l V3 # |cj ( j(&3) - ^(-^i)] (and permuta- 
tions) does not hold everywhere. 
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the singularity in the integrand in (4.10) is integrable, i.e., Bq is bounded. As a consequence we 
can perform the limit s — > + inside the integral. 



Proposition 4.6. The leading contribution Bq of the bispectrum of the fluctuations of the Bardeen 
potential *P has the form 

B (r,k u k 2 ,k 3 ) = -—- - 2 , 

1 2 3 

where Sq is a function ofk\T, k 2 T, and k^T only and k, = 

Proof. Analogously to proposition 4.3, we note that after a change of variables x — t p the 
integrand in (4.10) is a function of k\T, k 2 T, and k^T only. □ 

To finish our discussion about non-Gaussianities, we notice that, although the employed 
quantum field is a linear one, we obtained a three-point function for the Bardeen potential which 
is similar to the one obtained by Maldacena [25] who has quantized metric perturbations outside 
the linear approximation. 

5. Conclusions 

In this paper we have analyzed the influence of quantum matter fluctuations on curvature 
perturbations over de Sitter backgrounds. We used techniques proper of quantum field theory 
on curved spacetime to regularize the stress-energy tensor and to compute its fluctuations. In 
particular, we interpreted the curvature perturbations as the realization of a stochastic field. We 
then obtained the n-point functions of such a stochastic field as induced by the «-point functions 
of a quantum stress tensor by means of semiclassical Einstein equations. 

In order to keep contact with standard literature on inflation we allowed for curvature 
perturbations generated by simple metric perturbations related to the standard Bardeen potentials. 
We were motivated by the fact that, for classical scalar perturbation, standard Bardeen potentials 
O and *F are gauge invariant quantities and, whenever the perturbations are induced by matter 
described by a scalar field, the two are equal: <£ = *P. We thus assumed this point of view also in 
our semiclassical environment even if for stochastic or quantum fluctuations the latter assumption 
needs to be analyzed further because the corresponding gauge problem is more sever. 

We also noticed that, while the expectation value of the stress-energy tensor is characterized 
by renormalization ambiguities, this is no longer the case when fluctuations are considered. 
Hence the obtained results are independent on the particular regularization used to define the 
stress tensor. 

Although the considered model is certainly oversimplified to cover any real situation, it was 
possible to recover an almost-Harrison-Zel'dovich power spectrum for the correlation of the 
considered metric perturbation. Furthermore, the amplitude of such a power spectrum depends 
on the field mass which is a free parameter in our model and can be fixed independently of H. At 
the same time, since it does not depend on the Hubble parameter of the background metric, this 
indicates that it is not a special feature of de Sitter space. At least close to the initial singularity, 
the obtained result depends only on the form of the most singular part of the two-point function 
of the considered Bunch-Davis state. We thus argue that a similar feature is present in every 
Hadamard state and for backgrounds which are only asymptotically de Sitter in the past. 

Finally we notice that, since the stress-energy tensor is not linear in the field, its probability 
distribution cannot be of Gaussian nature. Thus we showed that non-Gaussianities arise naturally 
in this picture. 
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